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Two-phase  flow  regimes  are  reviewed  briefly.  A  physical  model,  which 
is  based  on  the  concept  that  droplet  diffusion  through  a  steam  boundary 
layer  is  the  limiting  mechanism  for  burnout  in  turbulent  flow,  is  de¬ 
scribed.  An  equation  is  derived  relating  burnout  to  other  parameters 
in  fog  flow.  With  simplifying  assumptions,  an  order-of-magnitude 
agreement  between  analysis  and  experimental  burnout  data  in  fog 
flow  is  demonstrated. 
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Nomenclature 

c  =  average  liquid  droplet  concentration, 
Itju/cu  ft 

c  =  liquid  droplet  concentration,  Ittn/cu  ft 
D  =t  equivalent  diameter,  ft 
f  =  friction  factor 

Gt  =  total  mass  flow  rate,  Ibju/hr-sq  ft 

h;fg  =*  latent  heat  of  vaporization,  Btu/lbm 

k_  =  film  coefficient  for  mass  transfer,  fph 
®  . 

Kq  =  kg/ug  =  dimensionless  mass  transfer  coeffi¬ 
cient 

m  =>  liquid  droplet  mass  current,  lbjj|/hr-sq  ft 
Re  =  Rejniolds  number 
s  3  hg/uf  =>  aLxial  slip  ratio 
Sc  =  Schmidt  number 

u  o  average  velocity  at  a  section,  fph 
X  =  average  steam  flow  quality  at  a  section, 
weight  fraction 

fgQ  =  burnout  heat  flux,  Btu/hr-sq  ft 
/I  =  viscosity,  Ibm/hr-ft 

Subscripts 

f  =  liquid  phase 
g  =  vapor  or  gas  phase 
w  «•  wall 

INTROKJCTION 

Burnout  presents  one  of  the  principal  limi¬ 
tations  In  the  design  of  liquid  or  wet  steam- 
cooled  nuclear  reactors,  rocket  nozzles,  and 
other  hlgh-speolflc-power  equipment.  Because  of 
this,  the  subject  of  burnout  has  received  consid¬ 
erable  attention  during  the  past  decade.  The 
present,  state-of-the-art  in  the  U,  S,  and  abroad 
Is  well  described  In  several  reviews  It 

Is  apparent  from  these  publications  that  today's 
knowledge  Is  not  sufficient  to  predict  burnout 
with  the  desired  degree  of  accuracy. 

Analytic  methods  for  predicting  burnout  have 
utilized  physical  models  which  have  generally  been 
based  on  visual  observations  of  pool  boiling. 

These  obsei^atlons  have  led  to  methods  which  are 
tied  to  bubble  dynamics  or  to  the  movement  of  li¬ 
quid  and  vapor  spikes  In  regular  geometric  cells. 
Since  they  are  based  on  pool-bolllng  observations. 
It  Is  not  surprising  that  these  analyses  correlate 
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Burnout  in  Turbulent  Flow  - 
A  Droplet  Diffusion  Model 

pool-bolllng  data  quite  well.  However,  they  fall 
for  flowing  systems. 

It  Is  suggested  In  the  present  paper  that  the 
controlling  mechanism  for  burnout  In  flowing  sys¬ 
tems,  particularly  at  high  Reynolds  numbers.  Is 
not  to  be  found  In  bubble  behavior  but  rather  In 
the  eddy  diffusion-limited  transport  of  liquid 
droplets  through  a  steam  boundary  layer  to  the 
heated  wall.  This  concept  appears  particularly 
attractive  for  fog  flow,  but  may  also  be  applicable 
to  flowing  low-quality  steam  and  even  subcooled 
liquid  systems. 

This  concept  grew  out  of  an  earlier  one  des¬ 
cribed  by  Clcchltt,  et  al  (J_)  .  They  postulated 
that  the  heat- transfer  area  Is  covered  by  a  thin 
liquid  layer  which  Is  continuously  fed  by  water 
drops  carried  by  the  steam  flow.  During  the  prep¬ 
aration  of  the  present  paper  our  attention  was 
called  to  theories  by  Vanderwater  (^) ,  Singh  (£) , 
and  Fauske  (10)  ,  who  similarly  proposed  that  bum- 
out  occurs  when  the  liquid  mass  flow  through  an 
annular  film  and  from  diffusing  droplets  Is  In¬ 
sufficient  to  support  the  heat  flux.  Reynolds  (11) 
proposes  a  correlating  factor  which  Is  based  on 
all  the  liquid  being  In  an  annular  film  at  the 
wall. 

In  this  paper,  two-phase  flow  regimes  are 
briefly  reviewed,  a  burnout  model  Is  described, 
basic  equations  for  btimout  in  fog  flow  are  de¬ 
rived  and  a  preliminary  comparison  is  made  between 
analysis  and  the  results  from  burnout  experiments 
in  the  fog- flow  regime, 

TWO-PHASE  PLOW  REGIMES 

Much  of  the  burnout  data,  needed  for  the  de¬ 
sign  of  nuclear  reactors,  has  been  obtained  ex¬ 
perimentally  with  high-pressure  water  In  the  qual¬ 
ity  region.  It  Is  Instructive  to  attempt  to  vis¬ 
ualize  what  the  appearance  of  the  two-phase  mix-  • 
ture  might  be  In  this  region. 

The  flow  of  two-phase  mixtures  has  been  ob¬ 
served  by  numerous  Investigators  under  isotheimal 
conditions  and  has  been  summarized  In  References 
{^)  ,  (1^)  and  (22)  •  The  observed  two-phase  flow 
regimes  are  well  described  by  the  names  frequent¬ 
ly  applied  to  them;  namely,  annular,  slug,  dis¬ 
persed,  bubble,  froth,  and  so  on.  There  Is  no 
uniform  agreement  on  terminology  or  even  on  the 
number  and  boundaries  of  regimes  that  do  occur. 
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Steam  Quality  (x),  weight  percent 

Fig.  1(a)  Steam-water  system  flow  pattern  chart. - 14.  7  psia;  - 800 

psiaj -  1500  psia.  (Based  on  reference  14) 


One  can  nevertheless  construct  a  flow-pattern 
chajpt  for  steam  and  liquid  water,  based,  for  In¬ 
stance,  on  Baker's  plot  ( 14)  ,  as  shown  In  Plg.l, 
Keeping  In  mind  that  the  boundaries  between  re¬ 
gimes  shown  as  lines  In  Plg.l  are  really  not  sharp¬ 
ly  defined  and  that  heat  addition  In  small  chan¬ 
nels  will  move  the  boundaries  towards  lower  qual¬ 
ities  and  lower  mass- flow  rates,  one  can  still 


reach  certain  over-all  conclusions. 

First,  It  will  be  noticed  that  at  mass- flow 
rates  larger  than  500,000  Ib/hr-sq  ft,  which  are 
of  Interest  to  nuclear  reactors  and  other  high- 
velocity  systems,  and  for  which  most  of  the  ex¬ 
perimental  burnout  data  have  been  taken,  the  fluid 
Is  either  In  the  bubble-froth  or  In  the  dlspersed- 
fog  flow  regime.  Bubble- froth  flow,  as  shown  in 
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steam  Quality  (x),  weight  percent 

rig.  1(h)  Steam-water  system  flow  pattern  chart.  — - 14.  7  psia;  - 800 

psiai -  15,  000  psia.  (Based  on  reference  14) 


Plg.l,  occurs  at  steam  qualities  of  less  than  ap¬ 
proximately  10  per  cent  depending  on  pressure  and 
probably  on  heat  flux.  As  the  name  indicates.  In 
this  regime,  the  fluid  Is  essentially  a  liquid 
containing  vapor  bubbles.  At  steam  qualities  be¬ 
tween  10  and  100  per  cent  and  high  velocities, 
the  fluid  Is  essentially  a  fog;  l.e.,  steam  with 
a  dispersion  of  liquid  droplets. 


In  the  fog-flow  regime  most  of  the  liquid  Is 
carried  In  droplet  form  .by  the  steam,  although 
under  Isothermal  conditions  a  thin  liquid  film 
can  be  observed  at  the  wall.  Since  there  are  no 
thick  layers  or  slugs  of  liquid.  It  Is  difficult 
to  visualize  the  growth  and  movement  of  bubbles 
In  fog  flow,  and  one  Is  forced  to  abandon  for 
this  regime  burnout  models  which  are  based  on  bub- 
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ble  dynamics.  A  droplet  diffusion  model  Is  pro¬ 
posed  In  Its  stead  as  follows: 

DESCRIPTION  OP  BURNOUT  MODEL 

For  Plovf 

In  fog  flow,  liquid  droplets  are  more  or  less 
uniformly  dispersed  in  the  hlgh-veloolty  steam. 
Under  high  heat  flirxes,  the  liquid  which  Is  cool¬ 
ing  the  wall  evaporates  and  must  be  replenished 
by  droplets  from  the  main  stream  if  burnout  Is  to 
be  avoided.  It  Is  proposed  that  this  replenish¬ 
ing  is  accomplished  by  eddy  diffusion  of  droplets. 
The  amount  of  liquid  which  can  be  transported  to 
the  wall  Is  limited  by  this  turbulent  diffusion 
process,  and  one  may  postulate  that,  to  a  first 
approximation,  burnout  will  occur  when  the  heat 
flux  becomes  greater  than  that  needed  for  the 
complete  evaporation  of  all  liquid  droplets  which 
diffuse  to  the  wall.  At  any  rate,  this  heat  flux 
represents  an  uppei  limit  for  burnout. 

In  this  model,  the  vaporization  of  droplets 
at  the  heated  wall  provides  a  diffusion  sink  and 
thus  establishes  a  concentration  gradient  for  a 
net  current  of  droplets  to  flow  from  the  main 
stream  to  the  wall.  A  similar  diffusion  process 
has  been  Investigated  by  Alexander  and  Coldren 
(15) .  They  performed  a  series  of  experiments  to 
elucidate  the  mechanism  and  to  measure  typical 
rates  of  deposition,  on  the  walls  of  a  straight 
duct,  of  small  water  droplets  suspended  In  a  tur¬ 
bulent  air  stream.  In  their  experiments,  the  thin 
liquid  film  at  the  wall  acted  as  a  diffusion  sink, 
and  the  major  resistance  to  droplet  transfer  from 
the  air  stream  to  the  duct  wall  resided  In  a  rela¬ 
tively  thin  layer  of  air  adjacent  to  the  wall. 

Prom  experiments  like  these  and  further 
analyses,  It  should  be  possible  to  deduce  the 
droplet-diffusion  rates  as  functions  of  pertinent 
parameters.  By  assuming  that  these  diffusion 
rates  are  controlling,  functional  relations  be¬ 
tween  burnout  and  these  parameters  can  then  be 
established,  which  In  turn  would  permit  the  re¬ 
liable  prediction  of  burnout  heat  fluxes  for  con¬ 
ditions  of  engineering  Interest,  Some  Initial 
steps  for  doing  this  are  Indicated  later  In  sec¬ 
tions  dealing  with  the  development  of  basic  equa¬ 
tions  and  comparison  of  data.  It  should  be  noted 
that  no  assumptions  need  to  be  made  about  the 
evaporating  process  Itself,  since  liquid  diffusion 
rates  are  assumed  to  be  controlling. 

Although  knowledge  of  the  evaporating  process 
Is  not  required  for  the  development  of  the  equa¬ 
tions,  It  Is  Interesting  to  speculate  as  to  what 
this  process  may  be.  We  presently  believe  that 
droplets  reaching  the  wall  evaporate  there  with¬ 
out  forming  liquid  films  which  are  thick  enough 


to  support  bubble  growth.  The  resulting  vapor 
flows  away  from  the  wall  Into  the  main  stream, 
thus  opposing  the  movement  of  droplets  towards 
the  wall.  In  cases  of  so-called  "slow  burnouts" 
(l6)  a  considerable  fraction  of  droplets  may  never 
reach  the  wall  at  all  but  may  evaporate  In  the 
slightly  superheated  steam  boundary  layer. 

Other  Flows 

One  ma;-  speculate  that  the  dlff-jslon  of  drop¬ 
lets,  rather  than  the  growth  of  bubbles.  Is  also 
controlling  for  burnout  in  low  quality  and  even 
In  subcooled,  turbulent  flows.  The  basis  for 
this  hypothesis  Is  that  near  burnout  fluxes  there 
Is  a  layer  of  vapor  near  the  wall  and  the  "bound¬ 
ary"  layer  Is  In  fog  flow.  Again  no  bubbles  can 
really  grow  at  the  wall,  Questions  related  to 
the  thickness  of  the  fog  layer,  the  effective 
liquid  concentration  in  the  main  stream  and  rad¬ 
ial  sprays  need  to  be  amswered  before  quantitative 
analyses  could  be  undertaken  to  predict  burnout 
fluxes  based  on  a  droplet-diffusion  model  for 
these  flows. 

DEVELOPMENT  OP  EQUATIONS  FOR  FOG  FLOW 

The  problem  of  mass  transfer  of  liquid  drop¬ 
lets  from  a  turbulent  stream  to  a  solid  boundary 
may  bi  viewed  mathematically  In  a  manner  analogous 
to  the  treatment  of  turbulent  momentum  exchange  and 
heat  transfer.  Implicit  In  such  a  mathematical 
description  Is  the  assumption  that  the  primary 
resistance  to  mass  diffusion  occurs  In  a  layer 
of  low  turbulence  adjacent  to  the  solid  boundary. 

One  may  define  a  film  coefficient  for  mass 
transfer  by  the  equation 

kg=  -  c^^)  (1) 

If  the  solid  boundary  under  consideration  Is 
a  heating  surface,  then  the  potential  heat-remov¬ 
al  capacity  of  the  system  can  be  related  to  the 
liquid-droplet  mass  current  by  a  heat  balance.  In 
line  with  the  previously  described  model.  Burn¬ 
out  may  be  assumed  to  occur  when  the  surface  heat 
flux  Is  equal  to  the  potential  heat-removal  capa¬ 
city  of  the  diffusing  liquid.  Combining  these 
concepts  we  obtain  the  basic  correlation  relating 
burnout  to  the  diffusion  properties  of  the  system, 
or 

?B0  -  *^ghfg(5  -  cj  (2) 

With  the  assumption  that  all  the  liquid 
reaching  the  wall  Is  evaporated  at  the  condition 
of  burnout,  the  wall  becomes  a  diffusion  sink  and 
c^,  >•  0,  Thus,  Equation  (2)  becomes 

^BO  “  *'g^fg°  (5) 
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The  average  liquid-droplet  concentratlori  In 
the  main  stream  can  he  approximated  by 

0  =  (1  -  x)G,./Uf  (4) 

This  assumes  that  the  droplet  concentration  is 
essentially  uniform  over  most  of  the  cross  section 
and  that  the  concentration  gradient  for  diffusion 
exists  In  a  relatively  thin  film  near  the  wall. 
Combining  Equations  (3)  and  (4)  and  multiplying 
the  right-hand  side  of  the  resulting  equation  by 
u_/u-  we  obtain 


^BO  77®  "  ^)'^t’^fg 

Ug  U^ 

(5) 

?B0  “  "  ^)®t'^fg 

(6) 

This  Is  the  basic  equation  for  predicting 
burnout.  It  cannot  be  readily  used  until  more 
Information  Is  available  on  mass- transfer  coeffi¬ 
cients  and  slip.  As  will  be  shown  later,  the 
principal  uncertainties  In  estimating  mass  trans¬ 
fer  coefficients  lie  in  a  lack  of  knowledge  of 
eddy  slip.  Thus  two  types  of  slip  are  of  concern; 
axial  slip  and  eddy  slip. 

Axial  Slip 

In  any  flowing  system  in  which  liquid  drop¬ 
lets  are  Introduced  Into  a  hlgh-veloolty  vapor 
stream,  the  droplets  will  experience  acceleration, 
or  deceleration,  until  their  velocity  approaches 
that  of  the  carrier  fluid.  In  a  fog- flow  heat- 
transfer  system,  the  carrier  fluid  is  continuously 
accelerating.  Thus,  unless  the  acceleration  of 
the  carrier  fluid  Is  occurring  at  an  extremely 
slow  rate  the  droplets  will  always  tend  to  lag 
the  carrier  fluid;  l.e.,  axial  slip  will  persist. 
Axial  slip  has  been  measured  at  a  number  of  In¬ 
stallations,  but  the  available  data  are  too  frag¬ 
mentary  to  be  of  practical  use  for  present  pur¬ 
poses,.  Levy  (12.)  reviews  these  data  and  proposes 
a  model  for  the  prediction  of  slip.  For  use  In 
Equation  (6)  we  suggest  that  slip  Is  a  function 
cf  a  number  of  parameters  such  as 

s  -  s  (P,a,x,  L,D)  (7) 

dJj 

and  that  considerably  more  work  needs  to  be  done 
to  elucidate  this  function.  Preliminary  studies 
indicate  that  axial  slip  Is  near  unity  In  fog 
flow. 

Eddy  Slip 

The  eddy  motion  of  the  carrier  fluid  has 
been  postulated  as  the  mechanism  which  transports 
liquid  droplets  from  one  point  in  the  stream  to 
another.  The  droplet  transport  to  the  wall  could 


be  calculated  from  an  analogy  between  naco  ani 
momentum  transfer  In  turbu-ent  streams  (_!£ >  • 

This  analogy  assumes  that  the  diffusing  material 
can  follow  the  eddy  motion  of  the  carrier  fluid 
exactly;  that  is,  eddy  dlffuslvltles  for  mass 
and  momentum  transfer  are  equal.  However,  in  the 
eddy  diffusion  of  liquid  droplets,  the  greater 
inertia  of  the  droplets  prevents  them  from  follow¬ 
ing  the  carrier  fluid  motion  exactly.  This  "eddy 
slip"  results  in  lower  values  for  mass  eddy  ’If- 
fuslvltles  and,  therefore,  mass  transfer  coeffi¬ 
cients,  than  would  be  obtained  from  the  analogy. 

An  estimate  of  this  effect  was  made  by  Long- 
well  and  Weiss  (_20)  for  the  eddy  slip  of  a  45  ^ 
diameter  kerosene  droplet  in  atmospheric  air  flow¬ 
ing  through  a  6-ln,  duct  at  a  velocity  of  3CG  fps. 
For  these  conditions,  they  found  the  ratio  of  the 
eddy  dlffuslvltles  of  mass- to-momentum  transfer 
to  be  0.35. 

COMPARISON  OP  ANALYSIS  AND  EXPE-RIMENTS  IN  KG  FLOW 


Although  we  do  not  kr.ow  the  exact  values  for 
axial  and  eddy  slips,  some  preliminary  oc.mfarlscn 
of  analysis  and  experiments  can  be  .made  tc  indicate 
that  the  proposed  model  does  show  promise. 

For  this  purpose,  one  may  make  the  assumptlcn 
that  the  axial  slip  ratio  is  equal  to  one  and 
calculate  mass  transfer  coefficients  from  burncut 
data  using  Equation  (6) 


-  r  '^BO 

"  =  1(1  -  xlGthfg^ 

Following  the  analogy  between 
turn  transfer,  Lin,  et  al  (22)  have 


(S) 

s  =  1. 

mass  and  momen- 
showr.  that 


Kq  =  il£  =  P(Re.  So)  (9) 

and  for  oases  in  which  the  Schmidt  number  is  near 
unity,  the  analysis  reduces  to 


Kq  =  =  I  =  P(Re)  (1C) 

Ug  2 

Thus  one  may  plot  mass  transfer  coefficients 
and  friction  factors  against  Reynolds  number  to 
test  the  validity  of  the  postulated  burnout  model 
by  observing  the  spread  of  data,  keeping  in  mind 
that  the  spread  may  be  considerable  because  of  the 
oversimplifying  assumptions  of  axial  slip  ratios, 
eddy  slip  ratios  and  Schmidt  number  equal  to  1. 

Fig, 2  serves  this  purpose.  It  shows  mass 
transfer  coefficients,  as  calculated  from  Equation 
(8) ,  plotted  against  Reynolds  number  for  all  ex¬ 
perimental  burnout  data  which  have  been  available 
to  us  In  the  steam  quality  range  of  2C  to  100  per 
cent  and  pressure  range  of  500  to  1500  psla.  The 
Reynolds  numbers  have  been  formed  by  the  arbitrary 
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Fig.  2(a)  Dimensionless  mass  transfer  coefficient  versus  Reynolds  number 
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use  of  saturated  steam  values  for  viscosities  and 
total  steam  and  liquid  flows  for  calculation  mass 
velocities. 

Also  shown  in  Pig. 2  are  the  droplet  mass 
transfer  coefficients  measured  directly  by  Alex¬ 
ander  and  Coldren  and  the  friction  factor  line, 

V2  =  0,023/Re®*^.  In  view  of  the  oversimplifying 
assumptions,  the  order  of  magnitude  agreement 
between  mass  transfer  coefficients  and  f/2  is  re¬ 
markable  . 

As  indicated  in  Pig. 2,  some  of  the  burnout 
data  had  been  taken  in  the  annular  and  slug  flow 
regimes.  Their  spread  is  large  because  the  sim¬ 
plifying  assumptions  which  are  more  nearly  applic¬ 
able  to  fog  flow  are  not  valid  in  these  flow  re¬ 
gimes. 

The  spread  of  the  fog-flow  burnout  data  is 
also  considerable,  but  can  be  explained  qualita¬ 
tively  by  slip  ratio  considerations.  It  is  hoped 
that  further  work  on  gaining  insight  into  turbu¬ 
lence  levels  of  two-phase  mixtures,  droplet  size 
distributions  and  droplet  diffusion  processes 
will  provide  functions  for  mass  transfer  coeffi¬ 
cients  and  slip  ratios  which  can  then  be  used  to 
predict  burnout  heat  fluxes  with  ♦^he  proposed  mo¬ 
del  . 

CONCLUSIONS 

1  It  is  proposed  that  droplet  diffusion 
through  a  steam  boundary  layer  towards  the  heated 
wall  is  the  controlling  mechanism  for  burnout  in 
fog  flow. 

2  It  is  suggested  that  the  same  mechanism 
nay  be  controlling  for  burnout  in  high-velocity, 
lew-quality  and  even  suboooled  flow  systems. 

3  An  equation,  relating  burnout  in  fog  flow 
to  variables  describing  the  droplet  diffusion  pro¬ 
cess,  has  been  developed. 

4  The  usefulness  of  this  equation  is  pres¬ 
ently  limited  because  data  for  certain  quantities 
such  as  axial  slip  and  droplet  mass  transfer  co¬ 
efficients  are  not  available, 

5  Applying  the  equation,  with  the  axial  slip 
ratio  equal  to  one*,  to  experimental  burnout  data 
in  fog  flow,  results  in  an  order-of-magnltude 
agreement  between  calculated  mass  transfer  coeffi¬ 
cients  and  single-phase  friction  factors. 

6  It  is  recommended  that  future  investiga¬ 
tions  of  burnout  in  turbulent  flows  concentrate 
on  gaining  insight  into  processes  related  to  the 
diffusion  of  droplets  and  their  evaporation  at 
heated  walls, 
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